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Abstract 

Using calculations from first principles we have investigated the lattice thermal conductivity 
of ideal mono- and bi-layer graphene sheets. Our results demonstrate that the intrinsic thermal 
conductivity of both mono- and bi-layer graphene is around 2200 Wm^^K^^ at 300 K, a value 
close to the one observed theoretically and experimentally in graphite basal plane, and at higher 
temperatures it decreases with the expected dependence. The little variation between mono- 
and bi-layer thermal conductivity suggests that increasing the number of layers does not affect 
significantly the in-plane thermal properties of these systems. 
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Although carbon-based low dimensional materials have been subject of intense investiga- 
tions for decades, the recent synthesis of graphene, a single 2-dimensional sheet of hexagonal 
carbon in the honeycomb lattice, has opened new avenues of research in the search of novel, 
alternative materials for microelectronics.^ The peculiar electronic structure of graphene^ 
implies that electrons have a negligible effective mass around the K point of the Brillouin 
zone (Dirac point) and thus a high electrical conductivity. However, electric currents are 
inherently coupled with Joule heating so it is imperative to investigate the thermal property 
of these systems to insure proper operations of an eventual device. Usually a good electric 
conductor is a good thermal conductor. However, this expectation is broken in the case of 
diamond, another kind of carbon-base crystal, which exhibit a high thermal conductivity re- 
gardless of its large band gap and low electric conductivity. Graphite, the closest relative of 
graphene in the family of carbon-based materials, displays high thermal conductivity along 
the basal plane and, not surprisingly, two orders of magnitude smaller thermal conductivity 
along the c-crystallographic axis.- 

There is a relatively large scattering in the experimental data for the thermal conduc- 
tivity of highly oriented graphite and graphene systems. Ref. O reports values of thermal 
conductivity in the range 1660 — 1880 Wm^^K^^. Balandin et air- measured values around 
4840 — 5300 Wm~^K~^ for a suspended graphene mono-layer. More recently, Nika et alA 
reported their theoretical investigation on the thermal conductivity of mono-layer graphene 
and publish values in a range between 3000 — 6500 Wm~^K~^, depending on the choice of 
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Griineisen parameters, 7, which had values ranging from 0.8 to 2.0 as deduced from Ref 
This particular study was based on the use of a valence force field method which resolves 
all interatomic forces into bond-stretching and bond bending modes.-^^ and the Griineisen 
parameters were assumed to be constant for each phonon branch.^ Comparable values of 
thermal conductivity have also been observed in Carbon Nanotubes (CNT): Kim et air- re- 
ported approximately 3000 Wm^^K"^ as a thermal conductivity of a multi-wall CNT with 
140 A diameter. Pop et alJ observed 2800 - 3900 Wm-^K-^ for single-wall CNT with 
diameters between 0.5 — 10 fim at room temperature. Both experiments were performed on 
a single suspended nanotube by heating it with electrical current. 

In order to resolve this apparent spread in the values of thermal conductivity and under- 
stand the interplay between the number of layers and the thermal response of these systems 
we have investigated the intrinsic lattice thermal conductivity of ideal mono- and bi- layer 
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graphene using phonon dispersion relations and lattice anharmonicity parameters calculated 
from first principles methods. Our results show that mono- and bi- layer graphene sheets 
have an intrinsic thermal conductivity superior to usual bulk solids except diamond which 
has 600 ~ 2000 Wm-^K-^ 4 

The lattice thermal conductivity k of a crystal at finite temperature T can be written as 

K = ^Kx (1) 

A 

«^A = $](v(q)-t)Mq)C^p.M (2) 
q 

where t is a unit vector in the direction of the thermal gradient VT, v(q) is the group 
velocity of the phonon modes, and r(q) is their life time.-i^>ii Here, the index A runs over 
the phonon modes and Cph{uj) is the contribution of phonon modes to the specific heat 
whose form is 

where fc^ is Boltzmann constant and h is Plank constant. In order to evaluate accurately 
the lattice thermal conductivity one needs to obtain the phonon group velocity v(q) from 
realistic phonon dispersion relations and to estimate the phonon life time from a careful 
consideration of the possible relaxation mechanisms. 

n 

As shown very elegantly in Ref 12|, one can calculate the graphite thermal conductivity 
along the basal plane considering only the LA and TA branches. The assumption is that the 
phonon dispersion of graphite has cylindrical shape and the ZA modes strongly interact only 
along the c-direction, mainly due to the large spacing and the weak bonding between the 
layers. Since the group velocities of the optical modes are considerably smaller than those 
of the acoustic modes, optical branches can be disregarded. Therefore, we can consider only 
TA and LA modes in the range above the temperature which corresponds to the highest 
ZA frequency and we can replace the first Brillouin zone (FBZ) with a circular cylinder 
and define the average sound velocity (f ) for the two-dimensional phonon gas consisting of 
LA and TA modes as: 2/(t>)^ = 1/{vla)'^ + '^/{vta)'^- Recent first principles calculations 
of phonon dispersion in graphite^*^ and our own results validate this approach since the 
phonon dispersion curves along the T — A axis are almost flat, implying that the FBZ of 
graphite looks like a circular column. Graphene sheets inherently have a two dimensional 
nature so it seems reasonable to accept these assumptions for the high temperature range in 



which the role of ZA modes is not significant. According to our calculations, the maximum 
of ZA mode is at 535 cm^^ . This frequency corresponds to 123 K {ksT = huj), validating 
our assumption for any T larger than that. 

We evaluated the average phonon life time as limited by the anharmonicity of lattice 
vibrations since this is the most fundamental limiting factor which is not relying on the 
purity of crystal or boundary termination technology. The analytical expression for the 
relaxation rate from the anharmonic three-phonon processes are readily derived as: 

1 = 27^^^ (4) 

where 7a are Griineissen parameters, ks is Boltzmann constant, T is temperature, M is the 
mass of atoms and A runs over the phonon modes.-i^^i^ Um and v are Debye frequency and 
the average sound velocity. 

The phonon dispersion relations used in this study have been obtained from first principles 
calculations within Density Functional Theory and Density Functional Perturbation Theory. 
The predictive power of this approach has been demonstrated by the excellent agreement 
between theoretical predictions and experimental observations in recent studies on the lattice 
vibrational modes of graphite and graphene.— In this study, the calculation for phonon 
dispersion and Griineisen parameters for the mode anharmonicity were obtained by using 
the PWSCF package of the QUANTUM-ESPRESSO distribution.i^ 

Fig. shows the calculated phonon dispersions along the high symmetry lines. The phonon 
frequency values of high symmetry points such as P, M, and K can be compared with previous 
studies and the results for mono-layer graphene are summarized in Table [B The results for 
bi-layer graphene show a very similar behavior. As shown in the figure, another optical 
vibrational mode appears at P of bi-layer graphene phonon dispersion. This ZO' mode at 
78 cm~^ is produced by the displacement of the carbon atoms in opposite directions along 
the c crystalline axis (optical mode). In graphite the same mode has a frequency of 95 
cm~^.— The agreement of our data with the previous theoretical and experimental studies 
is clearly excellent. From these accurate phonon dispersion curves we calculated the group 
velocity v(q) and the Griineissen parameter 7a for each phonon mode A. The latter are 
defined as the negative logarithmic derivative of the frequency of the mode with respect to 
volume: 7a = —[a/2u!\{q)][duj\{q)/da]^'^ where a is the lattice constant. By calculating 
phonon dispersion relation with small deviation from the original lattice constant and using 
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the above definition, we computed the Griineissen parameters over the FBZ as shown in 
Fig. . In mono-layer graphene, the ZO and ZA modes have negative 7a, between -1.38 and 
-0.17 for ZO, and between -53 and -1.46, for ZA. For the other modes, the values show 
variations between 0.16 and 2.76. Positive values of 7 correspond to a decrease of the 
frequency of the modes as the lattice parameter increases: the atoms are less tightly bound 
to their equilibrium position. Inversely, negative values of 7 correspond to an increase of the 
frequencies upon the expansion of the lattice. Since ZA and ZO are out-of-plane transverse 
polarization modes, one can imagine that if a graphene sheet is expanded, the force to move 
the atoms in the vertical direction will increase. As for the bi-layer, the dispersion relations 
for the Griineissen parameters display an additional contribution corresponding to the ZO' 
mode. In addition, 7za and 7zo' of of the bi-layer show slightly different behavior at the F 
point: they are positive. Small values of q (around F) correspond to long wavelength phonons 
and a positive value of 7 means that the atoms are less tightly bound. The variation of 7 from 
negative to positive tells us that the properties of the ZA or ZO' modes' collective movement 
are different for long wavelength and short wave length: as the volume is increased, and with 
it the intra-layer distance, ujza,zo' decreases at F. Atoms in both layers loose some of their 
coherence in the long wavelength vibrational modes. There are two TA and LA modes so 
the mean value of 7tai, 7ta2, 7lai, and 7la2 in the FBZ was found to be 0.52, 0.53, 1.58, 
and 1.56, respectively. 

The normalized thermal conductivity k/t is reported in Fig. . The curves are obtained 
from the ab initio phonon dispersion relations and the two-dimensional phonon gas model for 
infinite mono- and bi-layer graphene sheets along (100) direction from 200 K to 500 K. For 
mono- and bi- layer, we observe a similar behavior over the calculated temperature range: 
both show almost the same normalized thermal conductivities. For instance, at T=300 K, 
k/t = 4.78 X lQ^^Wm~^s~^ for mono-layer and n/r = 4.83 x 10^^Wm~^s~^. In both cases 
the contribution from the LA mode is larger than the one from the TA mode due to the 
large group velocity of the LA mode. From the data in Fig. (b), it is evident that the 
magnitude of the TA and LA modes is decreased when another layer is added to the system; 
however, also the number of modes doubles so that the normalized conductivity is the same. 
Moreover, going from mono- to bi-layer graphene, the slope of the TA and LA branches do 
not change significantly. However, the thickness of the layer doubles, so the contribution 
from each branch becomes one half of the value of the mono- layer.— 
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The lattice thermal conductivity with phonon mean life time limited by lattice anhar- 
monicity effects is presented in Fig. for the same temperature range of Fig and along the 
(100) direction. Both mono- and bi- layer graphene display a similar ~ 1/T dependence at 
high temperature. The values of k at room temperature (300 K) differ only by 8 Wm^^K^^ 
. As previously mentioned, the contribution from the LA branch is larger than that from 
the TA branch in the case of normalized thermal conductivity. However, as shown in the 
graphs, the TA mode transfers a larger portion of thermal energy than the LA mode in the 
real thermal conductivity. This is due to the longer phonon life time of TA which is expected 
from the comparison of the Griineissen parameters. The LA mode has three times larger 
7. Overall, nine times longer phonon mean life time is expected for TA so that this longer 
life time compensate the smaller group velocity. Klemens et a/.» predicted 1900 Wm~^K~^ 
as the lattice thermal conductivity along the basal plane for the bulk graphite at 300 K. 
Our results substantially agree with this previous estimate. The small difference, around 300 
Wm~^K^^ at room temperature, that we find can be explained by two reasons. First, in Ref. 
y they assumed 7=2 for both modes, a choice that eventually decreases the thermal con- 
ductivity. Second, they assumed an average sound velocity in the FBZ which could increase 
the thermal conductivity since the actual group velocity of phonon modes decreases when q 
approaches the zone boundary. The Debye frequencies Um of two dimensional phonon gas in 
mono- and bi- layer graphene are 1265 cm~^ and 1243 cm~^ respectively and those almost 
correspond to the maximum frequencies of the LA mode at the zone boundary which sup- 
ports the assumption that the frequency below Um (TA and LA) contributes to the thermal 
conduction. Finally, in Fig. we present the angular dependence of the thermal conduction. 
6 in the figure is the angle between the direction of the thermal gradient and the (110) 
crystalline direction. The angular dependence was calculated for different thermal gradient 
directions where, however, we observe only a negligible variation (less than 0.1% at 300 K). 

When it comes to realistic finite graphene nanoribbon, other possible scattering mecha- 
nisms can arise, in particular scattering by impurities and by the discontinuous structure of 
the ribbon side terminations. To evaluate these effects one should both measure the density 
of substitutional impurity atoms and the mass difference between the impurity atoms and 
the carbon atoms and gain a more complete understanding of the geometry of the system. 
Both pieces of information at not yet readily available from the current status of research 
on these systems. Therefore, our values for the intrinsic lattice thermal conductivity of 
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mono- and bi- layer graphene can be regarded as a theoretical upper limit of the thermal 
conductivity of more realistic graphene nanoribbons. 
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distribution, and a plane wave expansion up to 55 Ry. 32x32x1 Monkhorst-Pack meshes were 
used with a 0.02 Ry Fermi-Dirac smearing in the electronic occupations. We used the theoret- 
ical lattice parameter a =2.46A that is in excellent agreement with the experimental value at 
300K. The same simulation parameters were used for both mono- and bi-layer graphene setting 
the intralayer distance to the theoretical minimum 3.36A. The bi-layer graphene geometry is 
rhombohedral. 

We assumed a layer thikness h = 3.5 A for mono- 

layep22 and h = 3.5 + 3.36 = 6.86 A for the 

bi-layer. 
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Figure Captions 



Figure l(a,b). (Color online) Phonon dispersion of (a) a mono-layer and (b) a bi- layer 
graphene sheet along the high symmetry line. 

Figure 2(a,b,c). (Color onhne) Griineisen parameter of (a) a mono-layer and (b) a bi- layer 
graphene sheets along the high symmetry lines, (c) Two dimensional surface of Griineisen 
parameter of a mono-layer graphene sheet over the entire FBZ. 

Figure 3.(a,b) (Color onhne) Normahzed thermal conductivity in (100) direction of (a) a 
mono- layer and (b) a bi-layer graphene sheet. 

Figure 4.(a,b) (Color online) Thermal conductivity in (100) direction of (a) a mono-layer 
and (b) a bi-layer graphene sheet. Inset of (a): the angular dependence of the thermal 
conductivity of a mono- layer graphene sheet at T = 300 K. 6* is the angle between the 
thermal gradient and the (110) direction. 
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TABLE I: Comparison of the phonon frequencies of graphene at F, M, and, K in cm from various 
studies. 

Tzo Tlo/to Mza Mta Mzo Mla Mlo Mto Kza Kzo Kta Kla Klo Kto 
Thisstudy^o 884 1560 473 627 641 1318 1360 1396 532 550 997 1210 1228 1327 
Theoretical" 881 1554 471 626 635 1328 1340 1390 535 535 997 1213 1213 1288 
Experimental 861^ 1590'' 465*' 630^* 670** 1290^= 1321^= 1389^^ 482-^ 588'^ 1184'' 1184-^ 1313'^ 

'^Reference t2- Graphene sheets 

^Reference 12. Graphite 

'^Reference t? ■ Graphene sheets 

'^Reference 1 18l. Graphite 
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Fig. l(a,b): Kong et al. 
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Fig. 2(a,b,c): Kong et al. 
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Fig. 3(a,b): Kong et al. 
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Fig. 4(a,b): Kong et al. 
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